A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

06 (02.19.1; 08.03.4; 08.06.2; 09.08.1; 09.10.1; 13.09.6) 



ASTRONOMY 

AND 
ASTROPHYSICS 



ISO Spectroscopy of the Young Bipolar Nebulae S106 IR and Cep A 

East* 



M.E. van den Ancker^ 2, A.G.G.M. Tielens^^^, and P.R. Wesselius^ 



o 
o 
o 

Oh 

<: 



> ■ 

m ■ 
^ ■ 
o ■ 
o 
o 

Oh. 

6 ■ 



Astronomical Institute "Anton Pannekoek", University of Amsterdam, Kruislaan 403, NL-1098 SJ Amsterdam, Tlie Netherlands 

Harvard-Smitiisonian Center for Astropiiysics, 60 Gai'den Street, MS 42, Cambridge, MA 02138, USA 

Kapteyn Astronomical Institute, Groningen University, P.O. Box 800, NL-9700 AV Groningen, The Netherlands 

SRON, PO. Box 800, NL-9700 AV Groningen, The Netherlands 

NASA Ames Research Center, MS 245-3, Moffett Field, CA 94035, USA 

; accepted jdate^ 



Abstract. We present the results of ISO SWS and LWS grat- 
ing scans towards the embedded Young Stellar Objects (YSOs) 
SI 06 IR and Cep A East. Emission from the pure rotational 
lines of H2 and the infrared fine structure Unes of [C ll], [O l], 
[S l], [Si II] and [Fe ll], as well as absorption bands due to H2O, 
CO and CO2 ice were detected toward Cep A. In SI 06 we de- 
tected emission lines of H2, CO, H I, and a large number of 
ionized species including Fe, O, N, C, Si, S, Ne and Ar. S106 
also shows many of the infrared PAH bands in emission. Ex- 
citation temperatures and molecular hydrogen masses were de- 
rived from the low-lying pure rotational levels of H2 and are 
500 and 730 K and 8 and 3 x 10"'' Mq for S106 and Cep A, 
respectively. Since both objects are expected to have several 
solar masses of H2 in their environment, we conclude that in 
both cases the bulk of the H2 is cooler than a few hundred 
Kelvins. Excitation temperatures and line ratios were compared 
with those predicted by theoretical models for PDRs and dis- 
sociative and non-dissociative shocks. The [S l] 25.2 /im/[Si ll] 
34.8 /Ltm ratio is a particularly useful shock versus FDR dis- 
criminant and we conclude that S106 IR is dominated by PDR 
emission while Cep A East has a large shock component. From 
an analysis of the ionic lines in SI 06 we conclude that the cen- 
tral star must have a temperature around 37,000 K, correspond- 
ing to a spectral type of 08. From its luminosity it is concluded 
that the driving source of Cep A must also be a massive early- 
type star. The absence of strong high-ionization ionic lines in 
its ISO spectrum shows that Cep A has not yet created a signifi- 
cant H II region and must be younger than S 106, illustrating the 
process of the clearing of the suiToundings of a massive young 
star. 

Key words: Shock waves - Circumstellar matter - Stars: For- 
mation - H II regions - ISM: Jets and outflows - Infrared: Stars 
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1. Introduction 

The infrared emission-line spectrum of a Young Stellar Object 
(YSO) is dominated by the interaction of the central star with 
the remnants of the clouds from which it formed. The intense 
UV radiation generated by accretion as well as by the central 
star itself causes dissociation of molecular material close to the 
YSO and ionizes much of the atomic material, giving rise to 
typical nebular and recombination lines. Furthermore, the inter- 
action of the UV field with the remnants of the star's natal cloud 
will produce a photodissociation region (PDR; see HoUenbach 
& Tielens 1999 for a comprehensive review). In a PDR, heating 
of the gas occurs by collisions with electrons, photoelectrically 
ejected from grain surfaces. The gas in the surface regions of 
PDRs reaches temperatures of typically 500 K (e.g. Tielens & 
Hollenbach 1985), producing a distinctive infrared spectrum 
due to collisionally excited low-lying levels of many molecular 
and neutral species. 

Another phenomenon associated with YSOs are strong stel- 
lar winds, often collimated into a bipolar outflow. These will 
cause a shock as they drive a supersonic wave into the sur- 
rounding molecular cloud, heating it in the process. Shocks 
are usually divided into two distinct categories: J- or Jump- 
shocks, and C- or Continuous-shocks. In J-shocks viscous heat- 
ing of the neutrals occurs in a thin shock front in which radia- 
tive cooling is insignificant, while the post-shock gas is heated 
to several times 10^ degrees, dissociating all molecular ma- 
terial (e.g. Hollenbach & McKee 1989). Cooling of the post- 
shock gas occurs through atomic fine-structure lines as well as 
through re-formation of molecules. In contrast, C-shocks are 
magnetized, non-dissociative shocks in which the physical con- 
ditions change more gradually from their pre- to post-shock 
values and cooling is mainly through radiation from molecu- 
lar material (e.g. Kaufman & Neufeld 1996). If the tempera- 
tures in a C-shock become sufficiently high to start to dissoci- 
ate molecules, the cooling through the molecular lines dimin- 
ishes, and the shock temperature increases until it turns into a 
J-shock. Shocks with a shock velocity larger than 40 km s^^ 
are usually J-shocks, whereas slower shocks are usually of C- 
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Fig.l. SWS AOT SOI (solid rectangles), AOT S02 (dashed 
rectangles), AOT S06 (dotted rectangles) and LWS (solid cir- 
cle; beam average FWHM) aperture positions for our measure- 
ments of SI 06 superimposed on a K'-band image of the re- 
gion. The rectangles indicate the apertures (in increasing size) 
for SWS bands 1A-2C (2.4-12.0 /xm), 3A-3D (12.0-27.5 /im), 
3E (27.5-29.5 /xm) and 4 (29.5^0.5 /xm). 



type (Chernoff et al. 1982). In this paper we will illustrate the 
power of infrared spectroscopy to study the above mentioned 
phenomena using two young bipolar nebulae, S 106 and Cep A. 

The bipolar nebula SI 06 is one of the most studied Hii 
regions in our galaxy. The exciting source of the region is a 
very young massive stellar object, known as either IRS4 in the 
terminology of Gehrz et al. (1982) or IRS3 following Pipher 
et al. (1976). Around this central source, Hodapp & Rayner 
(1991) found a cluster of about 160 stars, embedded in the 
molecular cloud surrounding S106. Observed radio emission 
and Hi recombination lines from the region have been sug- 
gested to arise in a strong (M w 2 x 10"^ yr^^), fast (w 
200 km s~^), stellar wind, driving a shock into the surround- 
ing extended molecular cloud (Hippelein & Miinch 1981). A 
dark lane, largely devoid of any optical or infrared emission, 
separates the two lobes of the H ii region. This dark lane has 
been quoted many times as a prime example of a large (30") 
circumstellar disk, consisting of cool gas and dust, seen nearly 
edge-on (Bieging 1984; Mezger et al. 1987). However, Loushin 
et al. (1990) showed this structure to be an expanding ring 
of molecular material rather than a protoplanetary disk. Near- 
infrared ro- vibrational emission of H2 is seen to arise just out- 
side the H II region, suggesting an origin in a PDR (Hayashi et 
al. 1990). 




Fig. 2. Same as Fig. 2 for Cep A East. 



Cep A is another well-known site of recent star formation. 
It contains a luminous (2.5 x 10^ L©; Evans et al. 1981) far- 
infrared source, as well as several fainter infrared sources. The 
core of the Cep A region remains obscured at optical and in- 
frared wavelengths due to massive amounts of extinction, and 
is the source of an energetic, complex molecular outflow (Bally 
& Lane 1982). Radio-observations show that it contains two 
strings of about 13 ultracompact H 11 regions, arranged in a Y- 
shape (Hughes 1988). The source IRS 6a (Lenzen et al. 1984), 
located in the infrared nebula to the east of the core, is the dom- 
inant source of the eastern lobe at 20 \xva (Ellis et al. 1990), 
but polarization measurements show that the nebula is illumi- 
nated by the compact radio source HW-2 (Hughes & Wouter- 
loot 1984), located ^ 5" south of IRS 6a (Casement & McLean 
1996). The fact that it is not visible at 20 /im implies an ex- 
tremely large extinction towards HW-2. H2 emission occurs 
in both the eastern and western parts of the nebula (Doyon & 
Nadeau 1988). Hartigan et al. (1996) showed that the molec- 
ular emission to the east appears as a regular jet, whereas that 
to the west concentrates in shells, which they proposed to arise 
in wakes from bow shocks surrounding Herbig-Haro objects in 
the region. ISO observations of the western part of the region 
were discussed by Wright et al. (1996), who modelled the ob- 
served infrared fine structure lines as arising in a planar J-shock 
with shock velocity 70-80 km s~^. The same authors mod- 
elled the observed molecular hydrogen emission as arising in a 
combination of several C-shocks. This multitude of shocks in 
the region agrees with the explanation of Narayanan & Walker 
(1996), who reported the presence of multiple episodes of out- 
flow activity from the region. Recently, Goetz et al. (1998) pre- 
sented strong evidence for this scenario through near-infrared 
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H2 and [Fe 11] images of Cep A East showing several distinct 
regions containing shocks. 

In this paper we will present new Infrared Space Observa- 
tory (ISO; Kessler et al. 1996) data on infrared fine- structure 
and molecular emission lines from the central region of S106, 
centered on IRS4, and the eastern part of Cep A, centered on 
IRS 6a. These data give for the first time access to a large spec- 
tral range that was unaccessible from ground-based observa- 
tions and contain lines which are decisive probes of the excita- 
tion conditions. We will show that the observed emission-line 
spectrum can be well explained as arising in the combination 
of a PDR and a Hll region in the case of SI 06, whereas this 
emission is shock-excited in the case of Cep A. 

2. Observations and data 

ISO Short Wavelength (2.4^5 /um) Spectrometer (SWS; de 
Graauw et al. 1996) full grating scans ("AOT SOI") of 
SI 06 IR and Cep A East were obtained in ISO revolu- 
tions 335 (JD 2450373.740) and 843 (JD 2450880.050), re- 
spectively. Each of these observations took 2124 seconds 
of observing time. In addition to this, deeper SWS grat- 
ing scans on selected molecular and fine transition lines 
in the two objects ("AOT S02") were obtained in revolu- 
tions 134 (at JD 2450172.705; S106), 220 (JD 2450258.463; 
Cep A) and 566 (ID 2450603.706; Cep A). In revolutions 558 
(JD 2450596.103) and 580 (JD 2450617.838), scans at the full 
SWS grating resolution ("AOT S06") covering the wavelength 



ranges of 3.0-3.5, 4.0-6.8, 7.0-7.6 and 12. 1-16.5 /xm were ob- 
tained for S 106 as well. ISO Long Wavelength (43-197 /um) 
Spectrometer (LWS; Clegg et al. 1996) grating scans ("AOT 
LOl") of SI 06 IR and Cep A East were obtained in revolu- 
tions 134 (JD 2450172.726) and 566 (JD 2450603.681), re- 
spectively. Data were reduced in a standard fashion using cal- 
ibration files corresponding to OLP version 7.0, after which 
they were corrected for remaining fringing and glitches. To in- 
crease the S/N in the final spectra, the detectors were aligned 
and statistical outliers were removed, after which the spectra 
were rebiimed to a lower spectral resolution. 

The SWS full grating scans consist of twelve different grat- 
ing scans, each covering a small wavelength region, which 
were joined to form one single spectrum. Because of the vari- 
ation of the diffraction limit of the telescope with wavelength, 
different SWS bands use apertures of different sizes. This is 
illustrated in Figs. 1 and 2, where we show their relative posi- 
tion and size, overlaid on K-band (2.2 jim) images of S106 and 
Cep A East by Hodapp & Rayner (1991) and Hodapp (1994). 
For a source that is not point-like, one may see a discontinu- 
ity in the spectra at the wavelengths where a change in aper- 
ture occurs, which can indeed be seen in both the SI 06 and 
Cep A spectra. The relative discontinuities in SI 06 are close 
to the maximum possible values, indicating that the source is 
extended across the entire SWS aperture at the longer wave- 
lengths (33" X 20"). In Cep A the discontinuities are smaller, 
pointing to a smaller size of the far-infrared source. 
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Fig. 4. Same as Fig. 3 for Cep A East. 



Since both spectrometers on board ISO use entrance aper- 
tures that are smaller than the beam size, a wavelength depen- 
dent correction has to be applied to the standard flux calibration 
when observing extended sources. We determined these cor- 
rection factors by convolving the K-band images of SI 06 and 
Cep A East, with the beam profile and applied these to the data 
before further analysis. The maximum correction to the flux 
was 8% in the region around 40 /xm. If the sources are more ex- 
tended at the longer wavelengths or in specific lines than the K- 
band images we employed in estimating the diffraction losses, 
we will have underestimated the flux. The maximum error this 
could introduce in the flux calibration is f«15%. 

Figures 3 and 4 show the resulting SWS and LWS full grat- 
ing scans for S106 IR and Cep A East. Scanned fines and mea- 
sured line fluxes or upper Umits (total flux for line with peak 
flux 3a) are listed in Table 1. The errors listed in Table 1 in- 
clude the errors in the absolute flux calibration; errors in line 
ratios used in the remainder of this paper may be lower. Plots of 
all detected lines, rebinned to a resolution A/AA of 1500 with 
an oversampling factor of four (SWS), or to a resolution of 500 
(LWS), are given in Figs. 5 and 6. In the same figures we also 
show the line profiles expected for an unresolved point-source 
and an unresolved extended source filling the entire aperture. 
All detected lines in both SI 06 and Cep A are compatible with 
these profiles. For lines measured with the full SWS grating 
resolution (AOT S02 or S06), this indicates that they have a 
FWHM of less than w 150 km s-^ 



3. Solid-state features 

The SWS spectrum of S106 (Fig. 3) consists of a relatively 
smooth continuum, with numerous strong emission lines su- 
perimposed. The continuum appears to consist of two compo- 
nents, with the break between the two occurring around 10 /im. 
This is similar to the situation seen in other regions of recent 
star formation, where these two components are commonly as- 
cribed to emission due to dust close to the star (2-10 /im) and 
small dust grains in the wide circumstellar environment (10- 
200 /Ltm; e.g. Cesarsky et al. 1996). 

The familiar UIR bands at 3.3, 3.4, 6.2, 7.6, 7.8, 8.6, 11.3 
and 12.7 /xm, usually attributed to polycyclic aromatic hydro- 
carbons (PAHs), are present in emission and strong. As first 
demonstrated by ground-based spatially resolved CVF spec- 
troscopy of the 3.3 nm band by FelU et al. (1984), these fea- 
tures do not come from the central point-source IRS4, but origi- 
nate in the wide circumstellar environment. The band-strengths 
of the UIR features in our ISO spectra appear consistent with 
the ground-based measurements present in literature (Hefele & 
Holzle 1980; FeUi et al. 1984; GebaUe et al. 1985). 

An absorption band around 3.0 /im due to the O-H stretch 
mode of H2O ice, previously undetected in S106, is present in 
the SWS spectrum. Although partly filled in by emission from 
the 7.6, 7.8, 8.6 and 11.3 /im UIR bands, a 9.7 /im absorption 
feature due to the Si-O stretching mode in amorphous silicates 
also appears present in Fig. 3. The combination of the contam- 
ination by the UIR bands and the limited baseline can explain 
why this feature has not been detected from the ground. A weak 



Table 1. Observed and extinction-corrected line fluxes for S106 IR and Cep A East (in 10 W m ^). Also listed are predicted line fluxes from 
the Hn + PDR model (S106 IR) and the shock model (Cep A) discussed in the text. 
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Fig. 5. Detected lines in SI 06 IR. The velocities are heliocentric. The grey Unes show the instrumental profiles for a point source 
(narrower profile) and an extended source (wide profile) filling the entire aperture. 
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Fig. 5. (Continued) 

Table 2. Colimin densities of solid state absorption features towards SI 06 and Cep A. 
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dip around 18-19 /xm in our ISO data may be identified with 
the Si-O bending mode of amorphous siUcate, confirming the 
presence of siUcate absorption in S106. 

Since extinction in the continuum surrounding the 9.7 /itm 
feature is small compared to the extinction within this fea- 
ture, the extinction ^1^ at wavelength A across a non-saturated 
9.7 jum feature can simply be obtained from the relation A\ 
= — 2.5 log(J//o). Using an average interstellar extinction law 
which includes the silicate feature (Fluks et al. 1994), we can 
then convert these values of A\ to a visual extinction, result- 
ing in a value of Ay = 13™7 toward SI 06 IR. This value is in 
excellent agreement with that of 13™4 ± 2.7, derived towards 
stars in the S106 embedded cluster (Hodapp & Rayner 1991), 
or the value of 15™± 3™predicted from infrared Hi lines by 
AIonso-Kosta & Kwan (1989). It is also in agreement with the 
extinction of Ay = 12™ towards the northern lobe of the SI 06 
nebula (Felli et al. 1984). It is somewhat larger than that to- 
wards the southern lobe {Ay = 8™; Felli et al. 1984) and much 



smaller than that towards the central source {Ay =21™; Eiroa 
et al. 1979). 

For Cep A East the SWS spectrum (Fig. 4) is dominated by 
absorption bands from a variety of ices. The O-H bend and 
stretch modes of H2O at 3.0 and 6.0 /itm are strong, as are 
the 4.27 iim 1^0=0 stretch and the 15.3 ^m 0=C=0 bend- 
ing mode of CO2. The ^^C=0 stretch of CO at 4.67 /xm is 
clearly detected. The absorption band around 7 /itm, attributed 
to solid CH4 (Boogert et al. 1996, 1997) is present as well. 
The unidentified absorption feature around 6.8 /xm, also ob- 
served toward NGC 7538 IRS9 and RAFGL 7009S (Schutte et 
al. 1996; d'Hendecourt et al. 1996) is also present in Cep A. 
The 9.7 jim amorphous silicate feature consists of a very deep, 
saturated, absorption. From the non-saturated wings of this fea- 
ture we can again derive a visual extinction, resulting in a value 
of Ay = 270™ for Cep A East. This value is within the range of 
Ay = 75-1000™ extinction for the central source and nebula 
reported by Lenzen et al. (1984). RAHs appear absent in Cep A 
East. 
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Fig. 6. Same as Fig. 5 for Cep A East. 



From the integrated optical depth J T{v)dv of a non- 
saturated absorption feature we can compute a column density 
N using an intrinsic band strength Am. For H2O, CO, CO2 
and CH4 ices, values of Am were measured by Gerakines et al. 
(1995) and Boogert et al. (1997). For silicates. Am is taken 
from Tielens & AUamandola (1987), based on the lab mea- 
surements by Day (1979). Integrated optical depth and column 
density values are Usted in Table 2. The derived abundances of 
100:6: 15 for H20:CO:C02 in Cep A East are within the range 
of values observed in other Unes of sight (Whittet et al. 1996; 
d'Hendecourtetal. 1996). 

4. Hydrogen recombination lines 

A rich spectrum of H i recombination lines is present in the 
SWS data of S106 (Table 1). In contrast. Hi lines are ab- 
sent from the infrared spectra of Cep A East. It is not clear 
whether the Hi lines observed with the rather large (27" x 
14") SWS beam in S106 are dominated by the strong (M ~ 
10-5 Mq yr-i) stellar wind (Hippelein & Miinch 1981; Felh 
et al. 1984), or could be due to the extended H 11 region seen in 
Ha and Br7 (Bally et al. 1998; Hayashi et al. 1990; Greene & 
Rayner 1994; Maillard et al. 1999). 

A clue to what is the case could come from the observed 
line profiles. As can be seen from Fig. 5, the H I lines observed 
in SI 06 are wider than the instrumental line profile for a point 



source. Unfortunately, the effect of line broadening of a com- 
pact source with FWHM of a few hundred km s"^, as expected 
for the stellar wind (Felli et al. 1985; Drew et al. 1993), would 
be about the same as that of observing narrow recombination 
lines from an extended H 11 region. However, from the fact that 
our Bra line flux of 1.7 x 10"^^ W m~^ is identical to that ob- 
served by other authors in smaller apertures (Garden & Geballe 
1986; Persson et al. 1988; Drew et al. 1993), we conclude that 
at least for Bra the dominant component must be the compact 
wind. In view of the limited spectral resolution of our data, a 
full analysis of the wind structure is beyond the scope of the 
present paper. 

5. Molecular hydrogen emission 

Both S 106 and Cep A East are well-known sources of extended 
molecular hydrogen emission (Longmore et al. 1986; Hayashi 
et al. 1990; Bally & Lane 1982; Doyon & Nadeau 1988; Goetz 
et al. 1998). All ground-based H2 measurements in literature 
refer to the ro- vibrational transitions, with upper energy levels 
E{{3) /k above 5000 K. These do not probe the bulk of the H2, 
which is expected to be at much lower temperatures. However, 
in both S 106 and Cep A East we have detected many pure rota- 
tional lines of H2 (Table 1), which have much lower upper en- 
ergy levels and are therefore able to directly probe the physical 
conditions in the dominant chemical species. And since their 
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Fig. 7. H2 excitation diagrams for SI 06 (top) and Cep A East 
(bottom). ISO observations of pure rotational lines are indi- 
cated by the circles. Triangles and squares indicate ground- 
based measurements of ro- vibrational lines from literature. The 
dashed lines give the Boltzmann distribution fits to the low- 
lying pure rotational lines and all lines with upper level ener- 
gies above 5000 K. The solid Une shows the sum of both con- 
tributions. 



transition probabiUties are quite small, these lines are optically 
thin and the excitation temperature will be close to the kinetic 
temperature of the gas. 

From the H2 Une fluxes I{J) listed in Table 1 it is possible 
to calculate the apparent column densities of molecular hydro- 
gen in the upper J levels, averaged over the SWS beam, N{J), 
using N{J) = with A the wavelength, h Planck's 

constant and c the speed of light. The transition probabilities 
A were taken from Turner et al. (1977). Line fluxes were cor- 
rected for extinction using the average interstellar extinction 
law by Fluks et al. (1994). For S106 IR we adopted the value 
of Av = 13'."7 derived in Section 3. The fact that the 0-0 S(3) 
line of H2, with a wavelength near the center of the 9.7 /zm 
amorphous silicate feature, was detected in Cep A indicates 
that it cannot suffer from the same extinction as the continuum; 
the H2 emission must originate in a spatially separate region 
from the continuum. Therefore we adopt a much smaller value 
of Av = 2™, expected to be a reasonable value for an origin 
in either a shock or PDR, for the emission lines observed to- 
ward Cep A East. Since the extinction correction at these mid- 
infrared wavelengths is small, a mis-estimate of Ay by a few 
magnitudes will not affect our results. 

A useful representation of the H2 data is to plot the log of 
A'^( J) /g, the apparent column density for a given J upper level 
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Fig. 8. Theoretical relation between a) Trot(H2) and n for PDR 
models, b) Trot(H2) and Vg for J-shock models, and c) Trot(H2) 
and v,^ for C-shock models. 



divided by the statistical weight, versus the energy of the up- 
per level, taken from Dabrowski (1984). These plots are shown 
in Fig. 7. Also plotted in these figures are several measure- 
ments of H2 ro- vibrational lines in SI 06 and Cep A from liter- 
ature (Longmore et al. 1986; Tanaka et al. 1989; Bally & Lane 
1982). Although these measurements were taken at the same 
positions on the sky as ours, the beam sizes used were differ- 
ent and hence one may expect to see systematic differences in 
the derived specific intensities if the source of H2 radiation is 
more extended than the beam size. In fact only the S106 mea- 
surements by Longmore et al. (1986), taken with the smallest 
beam diameter, 12", differ systematically from the other mea- 
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surements, indicating that the extent of the H2 emitting region 
in SI 06 is probably somewhat larger than their beam size. This 
agrees well with H2 1-0 S(l) images of the region (Hayashi et 
al. 1990; Maillard et al. 1999). In Fig. 7 the Longmore et al. 
measurements of S106 were scaled to the Tanaka et al. data. 

The statistical weight g used in Fig. 7 is a combination of 
the rotational and nuclear spin components. We have assumed 
the high temperature equilibrium relative abundances of 3: 1 for 
the ortho and para forms of H2 (Burton etal. 1992).ForaBoltz- 
mann distribution, the points in the plot shown in Fig. 7 should 
form a nearly straight line whose slope is inversely proportional 
to the excitation temperature, while its intercept is a measure 
of the total column density of warm gas. Using the formula by 
Parmar et al. (1991) and the rotational constants by Dabrowski 
(1984) to fit our data points in the low-lying pure rotational 
levels to a Boltzmann distribution, we arrive at values of 500 K 
and 9 X lO^^ cm'^ and 730 K and 5 x lO^^ cm'^ for S106 
IR and Cep A East, respectively. Using the distances to SI 06 
and Cep A of 1200 and 690 pc (Rayner 1994; Mel'nikov et al. 
1995), this corresponds to total molecular hydrogen masses of 
0.008 and 0.003 Mq within the SWS beam. Using the alter- 
native distance estimate towards S106 of 600 pc (Staude et al. 
1982) would decrease the S106 H2 mass to 0.004 Mq. The fit- 
ted Boltzmann distributions are shown as the leftmost dashed 
lines in Fig. 7. The fact that the points for the ortho and para 
form of H2 lie on the same line proves that our assumption on 
their relative abundances is correct. 

As can be seen from Fig. 7, both the pure-rotational and 
ro- vibrational lines of H2 with upper level energies higher than 
5000 K deviate significantly from the leftmost dashed line. In 
both cases a Boltzmann distribution was fitted to these lines 
as well, shown as the rightmost dashed line in both figures. In 
the case of S106, the relative location of the 1-0 and 2-1 S(l)- 
lines in Fig. 7 are indicative of fluorescent excitation through 
UV pumping (Draine & Bertoldi 1996; Black & van Dishoeck 
1997). The fitted Boltzmann distribution to the higher energy 
levels has thus no physical meaning, but may still be useful to 
provide a simple parametrization of the relative population of 
the energy levels. In the case of Cep A, the lines may indicate 
the presence of a smaller column (3 x lO'^^ cm~^) of hot (a few 
1000 K) molecular hydrogen, in addition to the large column 
of warm gas. 

Employing predictions of H2 emission from PDR, J-shock 
and C-shock models by Burton et al. (1992), Hollenbach & Mc- 
Kee (1989) and Kaufman & Neufeld (1996), we determined the 
excitation temperature T^ot from the low-lying pure rotational 
levels from these models as a function of density n and either 
incident FUV flux G (in units of the average interstellar FUV 
field Go = 1.2 X 10"" erg cm-^ s^^ sr^^; Habing 1968) or 
shock velocity Vg in an identical way as was done for the ob- 
servations. The resulting relations between Trot and n or Vg 
are shown in Fig. 8. As can be seen from these plots, the PDR 
and J-shock models predict a fairly small (200-540 K) range 
of resulting excitation temperatures, whereas in the C-shocks 
this range is much larger (100-1500 K). Furthermore, we see 
that in the model predictions for shocks the resulting Trot does 
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Fig. 9. CO excitation diagrams for SI 06 (top) and Cep A East 
(bottom). The dashed lines give the fits of the Boltzmann dis- 
tribution to the data points. 



not depend much on density, whereas for PDRs it does not de- 
pend much on G, suggesting that once the mechanism of the 
H2 emission is estabHshed, it can be used to constrain Vg or n 
in a straightforward way. 

Comparing the excitation temperatures of 500 and 730 K 
for SI 06 and Cep A with those plotted in Fig. 8, we note that 
for SI 06 this falls well within the range of PDR- and C-shock 
model predictions, but are too high compared to the ones ex- 
pected from J-shocks. The observed bright PAH emission fea- 
tures (Sect. 3) and the atomic fine-structure line spectrum (Sect. 
6) point towards a PDR origin for the H2 emission in S 106. The 
higher temperature for Cep A can only be reproduced by the C- 
shock models. Therefore we tentatively conclude that a dense 
(> 10^ cm~^) PDR seems to be the best candidate to explain 
the observed H2 emission in S 106 IR and a slow (« 20 km s'^) 
non-dissociative shock can explain the observed warm column 
of H2 in Cep A East. Since the regions we are looking at proba- 
bly only fill part of the SWS beam, the absolute intensity of the 
H2 emission listed in Table 1 can also be reproduced by these 
same models by varying the beam filling factor. 

6. Carbon-monoxide emission lines 

In both S106 IR and Cep A East several ro- vibrational emis- 
sion lines due to gas-phase CO were detected in the long- 
wavelength part of the LWS spectra (Figs. 5 and 6). Similar 
to what was done for the H2 emission in the previous sec- 
tion, we constructed a CO excitation diagrams, using molec- 
ular data from Kirby-Docken & Liu (1978). They are shown 
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in Fig. 9. The temperature and column of CO resulting from 
the Boltzmann fit to this excitation diagram are 420 K and 5.6 
X 10^^ cm-2 and 340 K and 1.4 x 10^* cm-^ for S106 and 
Cep A East, respectively. Corresponding CO masses are 1 x 
10"^ Mq (S106 at 1200 pc; for the 600 pc distance estimate 
this mass estimate would decrease to 3 x 10~^ Mq) and 1 x 
10-3 Mq (Cep A East). 

The CO excitation temperature of 420 K for SI 06 is some- 
what lower than that found from the H2 Unes, in agreement 
with what is expected from a PDR, in which the CO emission 
arises in deeper embedded regions than the H2. We conclude 
that a PDR is the most likely candidate for the source of the 
gas-phase CO emission in S106. For Cep A East the CO exci- 
tation temperature of 340 K is much lower than that found from 
the H2 lines. This behaviour is hard to reconcile with an origin 
in the C-shock invoked to explain the H2 observations. How- 
ever, a J-shock, necessary to explain the ionic lines which we 
will discuss in the next section, may produce the observed CO 
emission while only emitting little warm H2. Since the LWS 
CO observations were made with a much larger beam than the 
SWS H2 observations, an alternative explanation may be that 
we are observing emission from two spatially distinct regions. 
The relatively large mass of warm CO as compared to that seen 
in H2 may also easily be explained in this scenario. 

The observed CO lines have critical densities of around 
10^ cm""^. Therefore the detection of these lines also implies 
densities of this order of magnitude or higher in the originat- 
ing region. For the shock seen in Cep A East this may not be 
unreasonable. However, assuming that these densities would 
exist in the entire SI 06 PDR would be implausible. If these CO 
lines do indeed arise in the large-scale environment of S 1 06, the 
PDR must therefore have a clumpy structure (e.g. Burton et al. 
1990). An alternative explanation of the presence of these lines, 
would be to identify the originating region with the surface of 
the extended disk-like structure surrounding SI 06 which could 
act hke a PDR. With the present data-set we cannot make a 
distinction between these two possibiUties. 

7. Atomic fine structure lines 

Important constraints on the physical conditions in the fine 
emitting region come from the observed fine structure lines. 
As a first step towards identifying the mechanism responsible 
for the observed emission we can look at the mere presence 
of certain lines. The observed lines with high ionization poten- 
tials in S106, such as fOiii], [Arm], [Neiii], [Nim], [Fein] 
and [Siv] can only originate in the Hll region surrounding 
S106. The PAH emission as well as the molecular lines ob- 
served towards SI 06 are indicative of the presence of a PDR as 
well. This PDR might contribute to the observed [Fe 11], [Ni 11], 
[Sin], [Oi] and [Cii] emission. We thus need to model the 
fine-structure lines in SI 06 IR as arising in the combination of 
an H II model and a PDR. 

We have used the photo-ionization code cloudy (version 
90.04; Ferland 1996) to generate model predictions for line 
strengths for the Hll region surrounding SI 06 IR, assuming 



a spherical geometry and constant hydrogen density through- 
out the region. We generated a grid of models in which the 
input spectrum, taken from Kurucz (1991) models for a stellar 
photosphere, and the electron density in the H II region were 
varied. The total luminosity of the model was fixed to the to- 
tal bolometric luminosity of S106 of 4.2 x 10^ Lq (computed 
from the spectral energy distribution using data from literature 
and assuming a distance of 1.2 kpc). The line ratios of the ob- 
served [O III] and [S iii] fine structure Unes are expected to de- 
pend mainly on density, and hardly on the temperature. Their 
behaviour in our H II region model together with the observed 
line ratios in S106 are shown in Fig. 10. From the fact that both 
line ratios agree on the thus obtained value for the electron den- 
sity in the Hii region, 1.3-2.5 x 10^ cm^ \ we conclude that 
our assumption of a constant density is reasonable. 

Several other fines ratios ([Arll]/[Arlll], [Ne ll]/[Ne III] 
and [S iii]/[S iv]) depend mainly on the effective temperature 
of the star and do not depend on density much. Their behaviour 
in the CLOUDY models for SI 06 is shown in Fig. 10 as well. 
From these plots, it can be seen that the source of ionizing pho- 
tons should have a temperature between 37,000 and 40,000 K, 
corresponding to a spectral type of 06-08, in good agreement 
with earlier determinations (Gehrz et al. 1982; Staude et al. 
1982). However, not all observed line ratios yield the same 
temperature for the central star Most likely this is due to the 
fact that the Kurucz models used for the input spectrum do not 
include the opacity shortward of the He 11 ionization limit due 
to the stellar wind, which is expected to be strong in the case 
of S106. The true temperature for the central source is there- 
fore expected to be around the lower range of temperatures 
deduced from the Kurucz synthetic photospheres, i.e. around 
37,000 K. For illustrative purposes, we also show the ratio of 
[O I] 63.2 /im/[0 III] 51.8 /im for our model H 11 region, show- 
ing that virtually all the atomic oxygen will be in the form of 
[O III] and that we can thus safely attribute all of the observed 
[O I] emission to the PDR. 

To be able to explain the observed intensity of 1.3 x 
10-2 erg s-^ cm-i sr-^ for the [Sill] line at 34.8 /xm in S106, 
the Tielens & HoUenbach (1985) PDR models require a den- 
sity higher than w 10^ cm-'^ and G > 10^ Gq. This regime can 
also reproduce the observed ratios of [Si 11] and the [Fe 11] and 
[O I] lines, although to reproduce both the exact strength of the 
[Si II] emission and the relative strength to that to the [Fe ll] 
Unes requires a [Fe ii] depletion of « 30% higher than the one 
assumed in the Tielens & HoUenbach models. An alternative 
explanation could be that the PDR doesn't fill the SWS aper- 
ture at 35 micrometer (33" x 20") completely, increasing the 
surface brightness of aU Unes. The predicted intensities of [S i] 
and [Fe i] are sufficiently low to be undetectable, in agreement 
with the observations. 

Towards Cep A East, only fine structure Unes that can be 
produced in shocks were observed. In contrast to both C- and 
J-shocks, PDRs do not produce significant quantities of [S l] 
emission (Tielens & HoUenbach 1985). Therefore this Une 
must be completely due to one or more shocks. C-shocks only 
contain trace fractions of ions and hence cannot explain the 
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Fig. 10. CLOUDY model predictions for various emission line ratios from the S106 H II region as function of electron density Ue 
and temperature of the ionizing star. The hatched regions between the dashed Unes show the error interval of the observed ratios 
in S106IR. The solid Unes show the computed ratios for log rig = 1.0-6.0, or T^g = 15,000-50,000 K. The heavy curves indicate 
the lines with log rig = 3.0 or Tcff = 40,000 K (best fit), (a) Behaviour of the [O III] 51.8 /im/[0 III] 88.4 /im Hne flux ratio, (b) 
The same for the [S III] 18.7 /xm/[S III] 33.5 /um ratio, (c) Same for the [Ar ll] 7.0 /xm/[Arlll] 9.0 /xm ratio, (d) Same for the 
[Ne II] 12.8 /in]/[Ne iii] 15.6 fjm ratio, (e) Same for the [S iii] 18.7 /m]/[S iv] 10.5 /xm ratio, (f) Same for the [O i] 63.2 ij,m/[0 iii] 
51.8 /xm ratio. 



observed [Fell], [Sin] and [Cii] emission. Hence the sim- 
plest hypothesis would be to try to explain the observed fine- 
structure lines in Cep A in terms of a J-shock model. To explain 
the [Si] surface brightness of 1.0 x 10"'^ erg s~^ cm~^ sr~^ 
requires a moderately dense (10^-10^ cm~^) pre-shock gas. If 
the shock does not fill the ISO beam completely, as is expected, 
this number will go up. The absence of [Fe i] emission does 
indicate that we are dealing with higher (10^-10^ cm~^) den- 
sities, so line ratios seem to provide more reliable constraints 



in this case. However, the aperture sizes for lines measured at 
different wavelengths are also different in some cases, mak- 
ing this method only reliably applicable to line ratios mea- 
sured in identical SWS or LWS apertures. Therefore the ra- 
tio of [Fe II] 26.0 (Um to [S i] 25.2 /zm might provide the most 
reliable indicator of physical conditions in the shock. To repro- 
duce the observed value of 0.26 with the Hollenbach & Mc- 
Kee (1989) J-shock models requires either moderately dense 
(10^ cm~^) gas with a low (30 km s~^) shock velocity or a 
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high (10® cm^ ^) density with a faster (« 60 km s"'^) shock. To 
have the predicted [Fe i] emission sufficiently low to explain 
our non-detections of this species, the moderately dense, slow 
J-shock model is required. This regime iilso reproduces the rel- 
ative [Si II] strength. 

As was discussed in the previous section, one or multiple 
J-shocks cannot reproduce the observed H2 emission (although 
a contribution to this can be expected). The success of the J- 
shock model in explaining the observed fine-structure Unes and 
the absence of PAH emission in Cep A East, excluding the 
possibility of a significant contribution from a PDR, leads us 
to pose that a combination of one or more J- and C-shocks 
must be responsible for the observed emission in Cep A. The 
presence of more than one type of shock in the region could 
be linked to the reported multiple episodes of outflow activity 
from the embedded source (Narayanan & Walker 1996). In the 
presence of both a J- and a C-shock, the [Si 11] and [Fe 11] emis- 
sion would originate completely in the J-shock, whereas both 
the C- and the J-shock would contribute to the observed [S l] 
and H2 spectra. With the J-shock parameters derived from the 
fine-structure lines, an additional C-shock component with a 
shock-velocity of about 20 km s~^ is required to reproduce the 
observed pure rotational H2 emission. The observed hot col- 
unm of H2 in Cep A may then be due to formation pumping 
(the effect that H2 gets re-formed with non-zero energy in the 
post-shock gas after being dissociated in the shock front) in the 
J-shock. To have the ro-vibrational lines of comparable inten- 
sity as the rotational lines (through collisional excitation) re- 
quires pre-shock densities of the order of 10® cm^"^. For both 
SI 06 IR and Cep A, predicted Unes strengths of the best fit 
models are also listed in Table 1. 

8. Discussion and conclusions 

In the previous sections we saw that the infrared emission-hne 
spectrum of SI 06 could be well explained by the presence of 
an H 11 region and a clumpy PDR. The average density in the 
H II region should be around 2 x 10'^ cm^'^. The density in the 
PDR, which was estimated to be 10^-10^ cm~^, is several or- 
ders of magnitude higher, in agreement with one would expect 
in a scenario in which a young massive star has photo-ionized 
its natal cloud in all but the densest clumps. 

In Section 7 we concluded that the incident FUV flux on the 
S106 PDR should be raflier high (> lO^^ Go). We also showed 
that the central source in S106 is of spectral type 08, with a 
total luminosity of 4.2 x 10^ Lq. From a Kurucz (1991) model 
for a stellar photosphere with T^s = 37,000 K and log g = 4.0, 
we compute that such a star emits a total FUV (6-13.6 eV) flux 
of 3.6 X 10^® erg s"! sr'^. To dilute this steflar FUV field to 
a value of 10^-10® Go, the PDR must be at a location 3-10 x 
10^ AU away from the central star, corresponding to a projected 
distance of 5-17". This projected distance is independent of the 
assumed distance towards S106. It is within the range allowed 
by the SWS entrance apertures and is compatible with an origin 
in the « 30" diameter region of H2 emission in the H2 1-0 
S(l) image of S106 by Hayashi et al. (1990). We conclude that 



a central OS star is sufficient to produce the FUV radiation field 
reaching the S 106 PDR. 

These results agree well with other determinations. Schnei- 
der et al. (2000) have spatially resolved the 157.7 /im [C ll] and 
63.2 jum [Oi] emission in S106 and concluded that there is a 
PDR region with density ~ 10^ cm~^ and strong (up to 8 x 
10^ Go) UV field in the immediate environment of SI 06 and 
a more extended dense (10^-10® Go) PDR region exposed to 
a lower intensity (300-500 Go) UV field. We are clearly dom- 
inated by their first component. This same component may be 
identified with the high-density (10^ cm"-^) fluorescent H2 re- 
gions identified by Hayashi et al. (1990). 

Both the presence of [Si] 25.2 /im emission and the H2 
spectrum point unambiguously to shocked gas as the origin of 
the infrared emission Unes in Cep A East. A combination of 
a slow (20 km s~^) non-dissociative shock and a faster (30- 
60 km s~^) dissociative shock are required to explain our ISO 
spectra. The density of the pre-shock gas in Cep A East is sim- 
ilar to that in the S106 PDR: 10^-10® cm"-''. Our J-shock com- 
ponent may be identified with the [Fe 11] 1 .64 /xm clumps ob- 
served by Goetz et al. (1998). The C-shocks could be due to 
their H2 knots which do not show ionic emission. 

It is interesting to compare the results derived here for the 
emission from the eastern lobe of Cep A with the results by 
Wright et al. (1996), obtained using the same instruments on 
board ISO, for the western part of the nebula. They derived a 
Trot of 700 ± 30 K for the H2 Unes with upper level energies 
up to 7000 K, and a temperature range up to 11 ,000 K for the 
higher upper level energies. They explained this H2 emission 
as arising from a combination of at least two C-shocks with 
different pre-shock density, shock velocity and covering factor. 
In addition to this, they also reported emission from [Nell], 
[Si] and [Sin], which they explained as arising in a planar 
J-shock with pre-shock density of lO'^-lO'' cm"'' and shock 
velocity 70-80 km s~^. Qualitatively, the detected lines in the 
eastern and western part of the nebula are well in agreement. 
The main difference between the observed lines seems to be 
that they are much more intense in the western part of Cep A, 
indicating that the densities there are a factor of 100 higher than 
those obtained for the eastern lobe, but the shock velocities are 
comparable in both parts of the nebula. The similarity between 
these two parts of Cep A are in agreement with a scenario in 
which the driving source of the molecular outflow, HW-2, went 
through multiple episodes of outflow activity, as suggested by 
Narayanan & Walker (1996). In this scenario, the J-shock com- 
ponent could be due to the most recent period of enhanced 
mass loss, whereas one or more C-shocks could be due to older 
generations of outflows. Alternatively, the stellar wind material 
could produce a (fast) J-shock while the surrounding molecular 
cloud might be swept up by a slower C-shock. 

As was already remarked by Staude & Elsasser (1993), the 
differences between the environments of these two massive em- 
bedded YSOs, SI 06 and Cep A, are remarkable. The mid- and 
far-infrared observations presented in this paper show that the 
difference between these two sources cannot be due to a dif- 
ferent orientation of their nebulae; it can only be a reflection 
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of their different evolutionary status, with Cep A being the 
younger of the two. In SI 06 the stellar wind and UV radiation 
of the exciting source have cleared and excited a sufficiently 
large region to create strong PDR emission, whereas in the case 
of Cep A, the central source is still heavily embedded and we 
only observe the interaction of the outflow with its surround- 
ings. In due time, Cep A will clear its surroundings, ionize 
the hydrogen, and evolve into a bipolar nebula quite similar 
to S106. 
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